High resolution Compton scattering as a Probe of the Fermi surface in the Iron-based 

superconductor LaOi_a;Fa.FeAs 
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We have carried out first principles all-electron calculations of the (OOl)-projected 2D electron 
momentum density and the directional Compton profiles along the [100], [001] and ]110] directions in 
the Fe-based superconductor LaOFeAs within the framework of the local density approximation. We 
identify Fermi surface features in the 2D electron momentum density and the directional Compton 
profiles, and discuss issues related to the observation of these features via Compton scattering 
experiments. 
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1. INTRODUCTION 

Since the discovery^ of superconductivity in a fam- 
ily of iron-based superconductors (pnictides), there have 
been a large number of studies of their electronic prop- 
erties that have revealed similarities between pnictides 
and cuprates. The theoretical prediction of a striped an- 
tiferromagnetic spin-density-wave (SDW) ground state^ 
was confirmed by neutron scattering'^. Superconductiv- 
ity is found in LaOFeAs with either hole doping^ or elec- 
tron dopingii^. First principles calculations^ find that 
the density of states (DOS) near the Fermi level {Ep) 
is predominantly due to Fe-c? orbitals. Owing to the ap- 
proximate S4 symmetry of the FeAs tetrahedra, these 
Fe-d orbitals split into lower lying e^ (da;2_y2 ,d3z2_r2) 
and higher lying i2g states {dxyAyzAzx)^- Theoretical 
calculations^ suggest that superconductivity may not 
be caused by electron-phonon coupling. Just as in the 
cuprates, the antiferromagnetic instability, which is sup- 
pressed by doping, is one candidate to explain unconven- 
tional superconductivity. Angle-resolved photoemission 
spectroscopy (ARPES) experiments have recently been 
carried out in BaFe2As2^'^°, a related pnictide. The su- 
perconducting gap in LaFeAs0i_a;F2: with x ~ 10% has 
been determined from the optical reflectance in the far- 
infrared region.'^ 

X-ray scattering spectroscopy in the deeply inelastic 
(Compton) regime provides a direct probe of the corre- 
lated many-body ground state in bulk materials while 
avoiding the surface sensitivity of ARPES. The use of 
modern synchrotron sourcesii makes it possible to inves- 
tigate complex materials via the measurement of direc- 
tional Compton profiles^. 

In this paper, we report first-principles computations 



of the 2D-projected electron momentum density (2D- 
EMD) and Compton profiles (CPs) in the iron-based su- 
perconductor LaOFeAs. We discuss Fermi surface (FS) 
images in the 2D-EMD and its anisotropy defined by sub- 
tracting a smooth isotropic function from the spectrum. 
Our analysis of the CPs reveals that FS features related 
to hole- as well as electron-pockets are more prominent 
in the CP when the momentum transfer vector lies along 
the [100] rather than the [110] direction. 



2. METHODS 

Our electronic structure calculatioirs are based oir the 
local density approximation (LDA) of density functional 
theory. An all-electron fully charge self-consistent semi- 
relativistic Korringa-Kohn-Rostoker (KKR) method 
is used"-. The compound LaOi_a;F2;FeAs has a 
simple tetragonal structure (space- group P4/nmm). 
We used the experimental lattice parametersii of 
LaOo.syFo.iaFeAs in which no spin-density-wave order 
was observed in neutron-scattering. A non-spinpolarized 
calculation was performed and the magnetic structure 
was neglected. Self-consistency was obtained for a;=0 
and the effects of doping x were treated within a rigid 
band model by shifting the Fermi energy to accommodate 
the proper number of electrons. The convergence of the 
crystal potential was approximately 10^^ Ry. The elec- 
tron momentum density (EMD) p{jpx,Py,Pz) was calcu- 
lated on a momentum mesh with step (Ap^^, Apj,, Ap^) = 
(l/16a, l/16a, l/16c)27r. The total number of points was 



14.58 X 10 within a sphere of radius 12.8 a.u. in mo- 
mentum space. The 2D-EMD p'^'^(px,Py) was calculated 
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FIG. 1: (Color online) (a) LDA band structure of LaOFeAs at fcz=0. (b) Theoretical 2D-EMD p^'^{px,Py) of LaOFeAs projected 
onto the (001) plane, normalized to p^'*(0,0). (c) 2D-LCW folded momentum density. Computed hole-like (green dots) and 
electron-like (blue dots) Fermi surfaces are marked. Note momentum is given in (c) in units of 27r/a, where a=7.6052 a.u. 



as 



P iPx,Py) = J p{px,Py,Pz)dpz 

while the Compton profile J{pz) is given by 
JiPz) = p{px,Py,Pz)dpxdpy. 



(1) 



(2) 



3. RESULTS AND DISCUSSIONS 

In Fig. [IJa), we show the LDA band structure of 
LaOFeAs. For x=Q, three bands cross the Ep around 
the r point, forming hole-like FSs [marked by green 
dots in (c)] while two bands cross Ep around M{-k,tt), 
forming electron-like FSs [marked by blue dots in (c)[. 
As electrons are added, the F centered FSs shrink and 
completely disappear around a;=0.13. The bands near 
Ep are dominated by the Fe d orbitals. The FeAs lay- 
ers are separated by insulating LaO layers, with the re- 
sult that the dispersion of these bands along F — Z is 
small and, apart from a small F-centered 3D hole pocket, 
the FSs are quasi two-dimensional. Based on the fully 
three-dimensional computations, we take advantage of 
this quasi two-dimensionality and investigate quantities 
in the kx — ky plane by integrating out the kz component. 

Figure[TJb) shows a map of the theoretical 2D-EMD^^ . 
This distribution can be described by an inverted bell 
shape with fourfold symmetry. The peak is at the zone 
center with tails extending over several unit cells. The 
dense contours around high symmetry points are signa- 
tures of the FS discontinuities. All these features are 
hidden behind the large inverted bell shaped signal. In 
order to investigate the Fermi surface topology in detail, 
we have employed both the 2D Lock-Crisp- West (LCW) 
foldingi^i^ and the 2D-EMD anisotropy. 



The 2D-LCW folding of the projected 2D-EMD 
p'^'^{px,Py) is defined by 



n'{kx, ky) — y ^ p^'^ikx + Gx, ky -f Gy), 



(3) 



Ga3 .G^i 



where n{kx,ky) gives the number of occupied states at 
the point [kx, ky) in the first Brillouin zone by summing 
over all projected reciprocal lattice vectors [Gx, Gy). The 
effect of the matrix element can be eliminated via the 
2D-LCW folding process of Eq. (3), which thus pro- 
vides a tool for focusing on the FS features. The the- 
oretical 2D-LCW folding shown in Fig. [TJc) has been 
smoothed using a Gaussian function with Ap = 0.17 a.u., 
which is typical of the resolutions available in high res- 
olution Compton scattering experiments. The positions 
and sizes of the FS pockets of the undoped parent com- 
pound LaOFeAs found by our KKR band calculation are 
shown as green dots for hole pockets and as blue dots for 
electron pockets. Before the application of the aforemen- 
tioned Gaussian broadening, n{kx, ky) shows a maximum 
nmax='29A and a minimum nmm=24.5. The difference 
nmax — nmin='i-d is Consistent with five bands crossing 
the Fermi level in the LDA calculation. Even after includ- 
ing experimental resolution, the FS features are still quite 
visible as seen in Fig. [IJc). The maximum of n{kx, ky) 
at M{tt,tt) is associated with the electron pockets; the 
minimum at F (0,0) is related to the hole pockets. 

Figure [2] shows the 2D-EMD anisotropy, found by 
subtracting a smooth isotropic function from the 2D- 
EMD. FS features show up as closely spaced contours 
around F(0,0) and M(7r,7r) points. The momentum den- 
sity around the T{orM) points in the higher zones is seen 
to be lower (or higher) than the average due to the pres- 
ence of hole pockets (or electron pockets). The zone-to- 
zone variation of intensity of these features can be under- 
stood as a matrix element effect associated with the sym- 
metry of the hole pockets at F and electron pockets at M. 
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FIG. 2: (Color online) Calculated 2D-EMD anisotropy in the 
parent compound LaOFeAs. Letters A and B label F points 
for odd and even fc-space sublattices of the hole-pockets as 
discussed in the text. The red lines are centered at M points. 
The rapid changes in the momentum density along these lines 
are FS signatures of the electron-pockets. 



For instance, the weak signal at the origin (0,0) can be 
understood since the bands crossing the Fermi level are 
predominantly d orbitals, whereas only an s orbital yields 
a significant contribution to the momentum density at 
the origin. Owing to interference effects, the FS features 
display a marked modulation from zone to zone. The 
Fe atoms in the unit cell are located at high symmetry 
positions, Fel (0,0,0) and Fe2 (0.5,0.5,0) (in units of lat- 
tice constants). The wavefunctions of these two Fe atoms 
show a constructive and destructive iirterference in mo- 
mentum space, which can be represented by the structure 
factor Sg = 1 + e-»'^(™+"), where G = {mx + ny)[27r/a] 
is a reciprocal lattice vector. Whereas Sg is largest when 
{rn + n) is even, Sg vanishes when (m + n) is odd. There- 
fore, the FS features show the alternating pattern seen 
in Fig. [51 For the FS hole-pockets centered at F, the 
FS features at B for odd (to -|- n) are much weaker than 
those at A for even (to + n) . Deviations from this rule 
are an indication of hybridization of the Fe orbitals with 
other orbitals. For the FS electron-pockets centered at 
M, a similar pattern is found. In Fig.[2]we show red lines 
crossing the B sites. The change of the momentum den- 
sity along the direction in which Ato = An is more rapid 
than chairges along a direction for which Ato = —An. 

The theoretical CPs along [100], [110] and [001[ are 
shown in Fig. [Ha). The CP along [001] is a smooth 
curve, since there are no Fermi breaks along this direc- 
tion. Thus, we can use this profile for highlighting FS 
features from the other profiles. 
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FIG. 3; (Color online) (a) Directional Compton profiles(CPs) 
J(p) along [100], [001] and [110[ for x=Q. (b)-(e) Differences 
AJ{p) between various pairs of CPs for undoped {x = 0) and 
doped (a; — 0.15) cases as indicated in the figures. Vertical 
arrows mark FS crossings with the electron pockets centered 
at M (blue) and the hole pockets at P (green) as discussed in 
the text. 

In Fig.[3)^b), we subtract the [001[ from the [100[ CP. 
The resulting periodic patterns occur at the P and M 
points and are associated with hole pockets (green ar- 
rows) and electron pockets (blue arrows) respectively. 
The CP has a dip within the hole pocket regions and a 
hump within the electron pocket regions. The FS breaks 
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are clearly visible and should be amenable to exploration 
via high resolution Compton scattering experiments with 
statistics high enough to numerically differentiate the dif- 
ference profiles. 

The same strategy is applied to the [110] CP in 
Fig. EI^c); however, the FS features are not as clear as 
in the [100] direction. The main reason is that the con- 
tributions of the hole- and electron-like FSs projected on 
[110[ overlap each other and tend to cancel out. The in- 
terference pattern acts to amplify this effect as follows. 
When the EMD is projected to form the CP, the projec- 
tions of A and B are distinct points along [110[. For hole 
pockets (centered at F), A (B) has strong (weak) signals 
associated with FS breaks. For electron pockets [centered 
at (tt, 7r)[, the red lines in the EMD are parallel to [110[ at 
the projection of A, while they are perpendicular to [110] 
at the projection of B. We notice that the FS signals are 
strong only along the red lines. Therefore, at the projec- 
tion of A (B) , the FS signals are strong (weak) for both 
hole and electron pockets. In simple terms, at A, strong 
hole pocket signals cancel strong electron pocket signals, 
while at B weak hole pocket signals cancel weak electron 
pocket signals. Thus [110] is not a suitable direction for 
studying the Fermi breaks. 

Next, we discuss how FS breaks disappear with elec- 
tron doping. The Fermi level for a;=0.15 shown in 
Fig. 1(a) at 0.075eV is obtained by a rigid band shift. 
At this doping level, all hole pockets around F are re- 
moved. As indicated by blue arrows in Fig. El^d) , the FS 
breaks associated with the FS electron pockets remain in 
the ]100] CP. Compared to Fig.[31Jb), the dips associated 
with the hole pockets have completely disappeared (green 
arrows). In Fig.g^e), we subtract the [100[ CP with 15% 
doping from the [100] CP with 0% doping. An interesting 
pattern of periodic maxima and minima appears around 



the F and M points which are identical with the posi- 
tions of the hole and electron pockets, respectively. This 
may prove the most promising method of detecting FS 
signatures. 



4. CONCLUSIONS 

In conclusion, we have identified FS signatures in the 
momentum density of LaOFeAs, finding alternating in- 
tensity patterns in the 2D-EMD due to the symmetry of 
the crystal. FS signatures for both hole- and electron 
pockets are shown to be relatively strong in the [100] CP 
in comparison to the [110[ CP. We thus conclude that 
the [100[ direction is the favorable one for studying FS 
signatures. Our analysis further indicates that a doping- 
dependent experimental study should be able to deter- 
mine at which doping level the hole-like FS's disappear. 
The present work sets a baseline for future experimental 
Compton scattering studies in the pnictides. 

5. ACKNOWLEDGMENTS 

This work is supported by the US Department of En- 
ergy, Office of Science, Basic Energy Sciences contract 
DE-FG02-07ER46352, and benefited from the alloca- 
tion of supercomputer time at NERSC and Northeast- 
ern University's Advanced Scientific Computation Cen- 
ter (ASCC). It was also sponsored by the Stichting Na- 
tionale Computer Faciliteiten (NCF) for the use of su- 
percomputer facilities, with financial support from the 
Nederlandse Organisatie voor Wetenschappelijk Onder- 
zoek (Netherlands Organization for Scientific Research). 



^ Y. Kamihara, T. Watanabe, M. Hirano, and H. Hosono, J. 

Am. Chem. Soc. 130, 3296 (2008). 
2 T. Yildirim, Phys. Rev. Lett. 101, 057010 (2008). 
^ C. de la Cruz, Q. Huang, J.W. Lynn, J. Li, W. Ratcliff 

II, J.L. Zarestky, H.A. Mook, G.F. Chen, J.L. Luo, N.L. 

Wang and P. Dai, Nature 453, 899 (2008). 

* H.-H. Wen, G. Mu, L. Fang, H. Yang and X. Zhu, Euro- 
phys. Lett. 82, 17009 (2008). 

5 G.F. Chen, Z. Li, G. Li, J. Zhou, D. Wu, J. Dong, W.Z. 

Hu, P. Zheng, Z.J. Chen, H.Q. Yuan, J. Singleton, J.L. 

Luo and N.L. Wang, Phys. Rev. Lett. 101, 057007 (2008). 
® L. Boeri, O.V. Dolgov, and A.A. Golubov, Phys. Rev. Lett. 

101, 026403 (2008). 

D.J. Singh and M.-H. Du, Phys. Rev. Lett. 100, 237003 
(2008). 

* C. Cao, P.J. Hirschfeld and H.-P. Cheng, Phys. Rev. B 77, 
220506 (2008). 

^ D.V. Evtushinsky, D.S. Inosov, V.B. Zabolotnyy, 
A.Koitzsch, M. Knupfer, B. Buchner, Viazovska, M.S., 
G.L. Sun, V. Hinkov, A.V. Boris, C.T. Lin, B. Keimer, 
A. Varykhalov, A.A. Kordyuk and S.V. Borisenko, Phys. 
Rev. B 79, 054517 (2009) 



H. Ding, P. Richard, K. Nakayama, K. Sugawara, T. 

Arakane, Y. Sekiba, A. Takayama, S. Souma, T. Sato, T. 

Takahashi, Z. Wang, X. Dai, Z. Fang, G.F. Chen, J.L. Luo 

and N. L. Wang, Europhysics Letters 83, 47001 (2008). 
^ M.J. Cooper, P.E. Mijnarends, N. Shiotani, N. Sakai and 

A. Bansil (editors), X-Ray Compton Scattering, Oxford 

University Press (2004). 
^ J. Laukkanen, K. Hamalainen, S. Manninen, A. Shukla, T. 

Takahashi, K. Yamada, B. Barbiellini, S. Kaprzyk and A. 

Bansil, J. Phys. Chem. Solids 62, 2249 (2001). (2004). 
^ A. Bansil, S. Kaprzyk, P.E. Mijnarends and J. Tobola, 

Phys. Rev. B 60, 13396 (1999). 
* Y. Qiu, M. Kofu, W. Bao, S.-H. Lee, Q. Huang, T. 

Yildirim, J.R.D. Copley, J.W. Lynn, T. Wu, G. Wu, and 

X.H. Chen, Phys. Rev. B 78, 052508 (2008). 
^ 1. Matsumoto, J. Kwiatkowska, F. Maniawski, M. Itou, 

H. Kawata, N. Shiotani, S. Kaprzyk, P.E. Mijnarends, B. 

Barbiellini and A. Bansil, Phys. Rev. B 64, 045121 (2001). 
^ D.G. Lock, V.H.C. Crisp, and R.N. West, J. Phys. F: Met. 

Phys. 3, 561 (1973). 



